Spectral Correlation and Response functions in Quantum Dots 
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We derive a general relation between correlators of density of states fluctuations and density 
response functions. It applies equally to quantum chaotic systems of pure symmetry (unitary, 
orthogonal, and symplectic) as well as to the crossover region between the universality classes. This 
relation is much more robust than Wigner-Dyson statistics; its validity extends to disordered metals 
with finite conductance and even to the Anderson insulators with large localization length. 
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One way to investigate a many-body system is to 
measure its response to an external perturbation. The 
response functions to which these measurements corre- 
spond can be presented through the two-particle Green 
function. In this Letter, we obtain and discuss generic 
properties of the response functions of quantum chaotic 
systems focusing particular attention on the behavior of 
disordered metallic grains (quantum dots). The exact 
analytical results which we obtain apply beyond the uni- 
versal regime described by the random matrix theory. 

Chaotic quantum dots and chaotic quantum systems 
in general are known to exhibit a large degree of uni- 
versality. Typically, this has been demonstrated by fo- 
cusing on the statistical properties of spectra. Perhaps 
the most widely studied statistical characteristic is the 
two-point correlator of density of states (DoS) fluctua- 
tions R 2 {u>) ||§]. Measuring the probability of finding an 
eigenstate at a distance uj in energy from a given state, 
it illustrates clearly such features as the repulsion of en- 
ergy levels, spectral rigidity, etc. In particular, for disor- 
dered metallic systems, an expression for R,2(ui) can be 
determined explicitly ||,|| and compared with the result 
from the random matrix theory. Recently the correla- 
tor was generalized to account for spectra that disperse 
as a function of some external tunable parameter [^|-^) . 
The parametric two-point correlator of DoS expresses the 
probability to find two states separated by an amount u> 
and x in energy and parameter space respectively. 

Rq,(ui) provides just one example of a two-particle aver- 
age Green function. Others which describe responses of 
the system to external perturbations involve the struc- 
ture of the eigenfunctions through matrix elements in 
addition to spectral properties. However, the striking 
feature to emerge from recent studies is the observation 



that the different response Green functions are, in fact, 
not independent, but connected by a set of differential re- 
lations (see Ref. fLOfl , and below). In this Letter, we will 
investigate the origin of these relations. Doing so, we will 
demonstrate that, for disordered quantum dots, these re- 
lations can be substantially generalized to encompass the 
non-universal regime of finite dimcnsionless conductance 
g. The exact analytical relations we derive apply beyond 
random matrix theory, where analytical expressions for 
different correlators are as yet unknown. Indeed the re- 
sults are valid even when the states are localized provided 
the localization length is sufficiently large. 

To make our discussion concrete, we focus on the prob- 
lem of a quantum dot subject to an external perturbation 
$ whose strength is controlled by a parameter X, 

H(X)^^-(p--A) 2 + U(r)+X$(r). (1) 

The random impurity potential U(r) is taken to be 
(^-correlated white-noise with (U(r)U(r')) = 8(r — 
r')/2irvT, where v denotes the average DoS and r is 
the relaxation time. A vector potential A is included 
in the unperturbed Hamiltonian and allows an interpola- 
tion between orthogonal (T-invariant) and unitary (bro- 
ken T- invariant) symmetry. The perturbation <&(r) is 
assumed to be local by comparison with the mean free 
path £, and to preserve the generic (spin-rotation/ time- 
reversal) symmetries of the Hamiltonian H(X = 0). 

The response of the metallic grain to the external per- 
turbation is embodied in the two-particle Green func- 
tion. The separation of fast and slowly varying modes 
(see below) allow the contributions to this average to be 
separated into three parts defined by three spatially de- 
pendent functions k, rip and nc (Fig. |]J): 
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[(Gi lXl (ri, r 2 ) G% X2 (r 3 , r 4 )> - G^(r 12 )G fl (r 34 )] 



2n 2 v 2 

= /(ri 2 )/(r 34 ) k(R 12l R u ; CI, X) + f(r 23 )f(r 41 ) n D (R 23 , R 14 ; CI, X) 
+f(r M )f{r ia ) n c (R 2 4,Ri 3 ;Cl,X), 



(2) 



1 



where G^(r 1} r 2 ) - {ri\{E - H(X) ± iO)' 1 ^) de- 
notes the Green function, O = E 2 — Ex, X = X2 — X\ , 

(r 4 + Tj)/2, and { 
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Tj\, and Rij 



rep- 



resents the ensemble average. The normalized Friedel 
function is defined by f(r) = Im Gq (r)/Im Gq (0) where 



G^ ,j4 (ri 2 ) = (G^'^(t"i, r 2 )) denotes the average Green 
function. For a d-dimensional system, f(r) can be pre- 
sented as a Bessel function, 

f(r) = T(d/2)(2/k F r) d / 2 -' J d/2 -x(k F r) e- r / 2e . 

We remark that our notation for the functions rip and nc 
is chosen to distinguish the Diffuson from the Cooperon- 
type contributions (Fig. 1). The decomposition implied 
by Eq. (||) relies on the rapid decay of f(r) on a scale 
comparable to the wavelength kp 1 <C £■ Since the typ- 
ical length scale of k and no,c l& t ne system size L, or 



the localization length £, Eq. (g) is applicable not only 
when the electronic states are extended, but even when 
they are localized, provided that L,£ ^> i. 

If the symmetry of the Hamiltonian itself is unchanged 
by the perturbation $, and H(X) belongs to one of the 
Dyson ensembles (pure symmetry cases), then nc = rip 
for the systems in which T-invariance is conserved (or- 
thogonal or symplectic class), while nc — for systems 
in which T-invariance is violated (unitary class). There- 
fore, according to Eq. (Q), for pure symmetries, all the 
statistical properties of the two-particle Green function 
are embodied in the two functions, k and n = njj. How- 
ever, remarkably they are not independent of each other. 
In particular, we will show that they are connected by 
the relation 



dr\dr 2 







dX 2 



$(ri) + $(r 2 ) \ 8_ 
2X ) dfl 



k(rx,r 2 ;n,X) 



d 2 
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drip(r) [n D {r,r;Q,X) +nc(r,r;fL,X)] 



(3) 



where 



<p(r) = /dridr 2 $(ri)/ 2 (ri2)$(r 2 ) 5(r 



ri + r 2 



)• (4) 



Eq. (|3j> represents the main result of the paper. Its range 
of validity (Q,t,X$t <C 1) encompasses a wide interval 
from the quantum chaotic regime on one hand to the lo- 
calized phase on the other. It implies strong constraints 
on the properties of disordered metals in regions where 
little is known about the explicit form of k and n. 

In confined or "zero-dimensional" cases where g — > 00 
and the coordinate dependence of k and n can be ne- 
glected (the random matrix limit), it makes sense to dis- 
cuss the volume averaged counterparts, 



k(Q,X) 
n DtC (n,X) 



drxdr 2 f fc(n,r 2 ;0,X) 
V 2 \n D ,c{rx,r 2 ;n,X) 



The former is related to the two-point correlator of DoS 
fluctuations through the identity R 2 (oj) = 1+Re[fc(w, 0)], 
while the latter represents the two types of density re- 
sponse functions. In the same limit an unfolding of the 
spectrum (ui = ft /A) by the average level spacing A, 
and a rescaling of the perturbation x 2 — C(0)X 2 , where 
G(0) = A- 2 ((dE a (X)/dX) 2 ) lifts the dependence of the 
correlators on the microscopic details of the Hamiltonian 
including the nature of the perturbation. This indepen- 
dence is a characteristic of systems which displays quan- 
tum chaos, and the universal correlation functions k(oj, x) 
and nc,D{u, x) can be applied quite generally. In partic- 
ular, studies of nn c have been used to investigate the 
response of chaotic systems, such as the dielectric re- 
sponse in complex periodic crystals, the quantum return 
probability in quantum dots W%i an( ^ the statistics of 



the oscillator strengths JTof| . In the universal limit, a di- 
rect evaluation of k{u),x) and np(u),x) established the 
zero-dimensional analogue of Eq. (||) [jl0| . 
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(5) 



Note that the linear term of $ in Eq. (||) is absent 
since the average drift of the spectrum is excluded by 
Tr$ = 0. Since Eq. (|J) relies on the statistical indepen- 
dence of wavefunction and spectra statistics, the differ- 
ence of Eq. (||) and Eq. (||) may imply the presence of the 
correlation between the wavefunction and the spectra. 

A similar result to Eq. (||) is found to apply beyond 
zero-dimension models when $ describes an irregular or 
random potential with some finite range. Since averaging 
over $ is expected, the linear term of $ in Eq. (^) van- 
ishes, and <p(r) is independent of r. This contrasts with 
the case of a Gaussian white noise random potential when 
the linear term does not vanish. In this case. Eq. (||) 
generates a relation between the return back probability 
P(t) defined by 



P(t) 



dQ e 



-iClt 



(Gl x (r,r)G^ +ntX (r,r)). (6) 



and DoS correlator k |12| 



VTp 



S(t,X = 0) = --^t z P(t), 



(7) 



where we introduce the form factor S(t, X) — 
A- 2 Je mt k(Sl,X)dn. 

Finally, before turning to the derivation of Eq. (|^) , let 
us examine what the formula implies for the relation be- 
tween the universal response functions in the crossover 
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region between the universality classes. If a T-invariant 
Hamiltonian is subject to a perturbation which violates 
T-invariance, the degrees of freedom which rely on the 
interference of time-reversed paths — the Cooper chan- 
nel — are suppressed while the Diffuson channel is un- 
affected This is reflected in the different behavior 
shown by the correlation functions njj and nc- 

The effect of a symmetry breaking perturbation can be 
demonstrated by generalizing the Hamiltonian in Eq. (|l|) 
to include two types of perturbation: 



H(X , X u ) =H + X $ + X u • $ u , 



(8) 



where $ G preserves T-invariance, whereas $ u = 
breaks it. Since there are, in principle, many different 
ways to violate T-invariance, we designate <fr u as a vector 



d 



± 



d 



1 d 



cos 9 



d{x 2 u ) 4 5(i 2 ) 



describing different possible perturbations and assume 
oc Sij. In general, the two-particle aver- 
age Green function (G% uXoi Xui G4,x o2 ,x u2 ) depends, 
not on each of the many parameters, but only on five: 
the energy difference, Xp = Xp 2 — Xpi (j3 = u,o), 
\X\ = \X u2 + X u i\/2 and an angle 8 specifying the "rel- 
ative direction," X X = \X\\X\ cos (9 @. 

In the universal or random matrix limit, the rescalings 

xl = C o (0; 0)A o 2 , x\ = C u {0; <x)X 2 u , x 2 = C u {0; co)A 2 
C (O;X) = -^((dE Q /dX p f), ((3 = o,u), 

lead to a generalization of Eq. (@) for the crossover inter- 
val between orthogonal and unitary ensembles: 



I 

Ax 2 



dcost 



1JP_ 
2~{h? 



(10) 



Since k(uj, x a , x, x u , cos 9) is known explicitly JIJ], this 
formula determines the response functions nc, ud within 
the crossover region. 

We now outline the main elements involved in the 
derivation together with that of the decomposition for- 
mula of Eq. (Q) . For brevity we will confine our discus- 
sion to the properties of T-invariant systems in which the 
symmetry is not violated by the perturbation. An exten- 
sion to the crossover regime can be achieved straightfor- 
wardly. Our approach and discussion is based on the field 
theoretic construction developed to describe the low en- 
ergy "excitations" or diffusion modes in disordered met- 
als, where statistical properties of two-point correlators 
are presented through an effective nonlinear cr-model in- 
volving 8x8 supermatrix fields Q [||. 

Following the conventions of Ref. the 8x8 super- 
matrix can be decomposed into 4x4 submatrices Q ab 
(a,b = 1,2) which reflect the advanced and retarded 
subspace. The spatially dependent correlations functions 
can then be expressed as 



2k(r u r 2 ;il,X) = 
2n D (r 1 ,r 2 ;n,X) -- 
2n c (rx,r 2 ;Q,X) = 



-l-<Q&(ri)Q§!(r a )) Q , (11a) 
-<^3(n)Qi(r 2 )> Q , (lib) 
-(Qll(ri)Q 2 A(r 2 )) (11c) 



where the integration is restricted to the manifold on 
which Q 2 = 1 3 and (• • -) Q = J DQ (■■■) e ~ F ^l The 
effective action F[Q] is given by 

F[Q] = Jdr STr {D(VQf + 2i[Q - A$(r)]QA} 

2 2 p 

-— j-^-A 2 J drip(r)STr[QAQA], (12) 

where D = v^r/d denotes the diffusion coefficient, and 
A = diag(l4, — 14) is the constant matrix that breaks 



the symmetry between the retarded and advanced Green 
functions. The explicit structure of Q corresponding to 
the three universality classes can be found in Ref. || . 

The decomposition formula of Eq. (||) is obtained by 
presenting the average two-particle Green function in 
terms of the supermatrix Green functions, 

(Glx>i,r 2 )Gf 2 , X2 (r 3 ,r 4 )) = {Qll{r x ,r 2 ) 0» (r 3> r i )) Q 
+ (Qllirx^) gg(r 3) r 2 ) + 0£(n, r 3 ) 0£(r 4 , r 2 )) Q , 

(13) 

where Q = [E -p 2 /2m - (O - A$)A/2 - iQ/2r}^ 1 . For 
small values of the parameters, the explicit dependence 
of Q on and X can be neglected. Taking account of 
the slow variation of Q(r), Q can be separated into fast 
and slowly varying spatial modes: 

Q(r u r 2 ) « G^(n 2 ) + i^f(r 12 )[Q(R 12 ) - A]. (14) 



Substitution of Eqs. (|13|jlJ) into Eq. (|ll|) leads to Eq. 
To obtain Eq. (|), we examine the derivative 



Afc^X) = Jdrd ri dr 2 (Ski + 



Ski), 



(15) 



where 

Sk x = ^ {$(r)STr{Q(r)A]Qll(r 1 )Q 22 (r 2 )) l 
v 2 v 2 X 



Sh> = — 



1G 



{^r)STr[(Q(r)A) 2 ]Qll( ri )Qll(r 2 )) 



Q 



Since Ski is given through $>(r)(5/5§(r))k(ri, r 2 ; £1, X), 
the first term can be recast in the form 



1 



dr x dr 2 $(n) 



d_ 

'dE x 



*(r a ) 



d 
dE 2 



k{r 1 ,r 2 ) 



dr 1 dr 2 [ \—k(n,r 2 ). (16) 
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To calculate (5fc 2 , we work with the symmetrized expres- 
sion for k and nu c m the fermion-boson subspace Jig], 



k 



1 



1 

32 



<STr[fc s Q n (ri)]STr[fc s Q 22 (r 2 )]) , (17) 



1 



(18) 



n D + n c = — <STr[fc s Q 12 (r 1 )fc s Q 21 (r 2 )]) Q , 

where k s = diag(l2, — 12), and utilize the identity, 

(STr[g 12 Q 21 (r)]STr[fc s g 11 (r 1 )]STr[fc s g 22 (r 2 )]) Q 

= (STr[Q 11 (r 1 )]STr[g 22 (r 2 )]STr[fc s Q 12 fc s Q 21 (r)]) Q . (19) 

The identity originates from the freedom to choose the 
fermion or boson field to calculate correlators, i.e., the 
fermion-boson supersymmetry of the tr-model. (This 
property was presented explicitly for the two-point cor- 
relator of Q- matrices in Ref. Jll|.) Note that since this 
fermion-boson symmetry is required to impose the cor- 
rect normalization, it is present even in the crossover in- 
terval between universality classes. Since STr[Q n ] = 
— STr[Q 22 ] by the saddle point condition on the Q- 
matrix, the right hand side of Eq. ( |l9| ) coincides with 
the second derivative of nu + with respect to 17. In 
this way, the second integral over <5fc 2 can be shown to be 
equal to the right hand side of Eq. (||) . 

We emphasize that, since the proof of Eq. (||) relies 
only on the supermatrix nonlinear a-model, it applies 
with the same generality. This implies that the identity 
encompasses the behaviors both of the mobility edge as 
well as in the localized phase. 

In conclusion, we have derived and discussed a rela- 
tion between spectral correlators and response functions 
in disordered metals which applies over a range which 
encompasses quantum chaos (random matrix limit), the 
metallic phase (large conductance) and even the Ander- 
son insulating regime with a large localization length. 
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FIG. 1. (a) Diagram representation of the two-particle 
Green function {G A (r 1 ,r 2 )G R (r a ,r 4 )). (b) DOS fluctua- 
tion-type contribution k. (c) Cooperon-type contribution nc- 
(d) Diffuson-type contribution hd- 
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